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with the precise data obtained at Saclay is achieved . 

PACS: 21.45+v, 24.70+s, 25.10+s 



Keywords: NUCLEAR REACTION d(p,p,p)n E=2.0 GeV, calculated 
vector analyzing power A y , tensor analyzing power A yy , polarization 
of the proton, covariant formalism, double-scattering, delta-excitation 



*e-mail: olegreb@mail.desy.de 



1 



Introduction 

The exclusive break-up of the polarized deuteron by the protons 
at a deuteron beam energy of 2 GeV was studied at Saclay [TJ [2] • 
The experiment was motivated by the idea to explore the impor- 
tant details of the deuteron wave function at the high internal 
momenta. As compared with the series of previous dp —> ppn 
experiments 0, HI 13 El E|, this experiment was designed to 
be sensitive to the ratio between the S- and D-components of 
the deuteron wave function brought by polarization observables. 
The sensitivity to the S/D ratio is maximum when the Impulse 
Approximation (IA) is valid. In this approximation one of the 
deuteron's nucleon is a spectator and another one interacts with 
the beam/target proton. As theory predicts, and it is confirmed 
experimentally, the IA is valid only when the internal momen- 
tum of the deuteron does not exceed « 200 MeV/c. An at- 
tempt of straightforward experimental testing of the deuteron 
wave function above this momentum could be justified only if 
other mechanisms beyond of the IA are well under the control. 
Therefore a good knowledge of other interaction mechanisms is 
of the great importance in the experiments of such kind. 

That is why the theory of the correction to the IA has a his- 
tory as long as the experiments initiated by this approximation. 
Chew and Goldberger [8J represented the scattering of the ele- 
mentary particles by complex nuclei as the multiple-scattering 
series T = J2kLiT^ k \ where k is the number of two-body NN 
interactions. The corresponds to the single-scattering and 
contains the IA term, corresponds to the double-scattering 
and so on. Analogous development follows from the Faddeev 
equations p. Everett (TO] was the first who scrutinized the 
double-scattering terms. He used the spin-dependent NN am- 
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plitude from the phase shift analysis (PSA) corrected so as to 
allow one nucleon to be off-shell and calculated the loop inte- 
grals of the double-scattering graphs, taking the NN amplitudes 
out of the integral at a fixed point. Golovin et al. jTT] used the 
similar method with the up-date NN amplitudes. 

Another approach was developed by Wallace . He took an 
advantage of the Glauber cancellation by the high-order terms in 
multiple-scattering series of those pieces of the double-scattering 
loop integral, which originate from off-shell states (T3]. This 
approach obtained an application by Punjabi et al. [Hj. 

The question about the role of the A-isobar in the (p,2p)- 
reactions often arises when discussing the discrepancies between 
the experiment data and predictions with the only NN rescat- 
terings. Yano [TI)j calculated the contribution of the A-isobar 
to the deuteron break-up reaction using Feynman diagram ap- 
proach but the corresponding amplitude was added incoherently 
to the other amplitudes. 

The analysis of the copious and precise data on the polar- 
ization parameters in exclusive dp — > ppn reaction, obtained at 
Saclay [2], requires to deal with the interference between the 
various amplitudes of the multiple scattering. In this paper we 
present such a model. It is based on the approach of Everett 
and includes the spin and isospin variables. The model takes into 
account the double-scattering contribution and also the graphs 
including the process 7r*d — >NN, the virtual pion being emitted 
by the beam/target proton. Since the ird ^NN reaction at these 
energies goes mainly through AN in the intermediate state, the 
graphs considered by Yano are also added in a coherently to the 
other contributions. Furthermore, the pion-nucleon scattering 
in other partial waves (S, P, D) is also considered. 
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When calculating the nucleon-nucleon rescattering contribu- 
tions to the deuteron break-up amplitude one needs to transform 
the NN matrices from nucleon-nucleon center of mass frames to 
the common laboratory frame. These transformations depend 
on the spin basis chosen for the one-nucleon states. Usually 
the canonical or helicity basis are used which are transformed 
with unitary two by two matrices depending on a nucleon mo- 
mentum. They are so called Wigner rotations and one needs 
to make the four different rotations corresponding to the two 
incoming and two outgoing scattering nucleons. The covariant 
basis fTBl [T7j transforming independently of the particle mo- 
mentum by the unimodular two by two matrices is free from 
this deficiency. An amplitudes in the covariant basis are usually 
referred to as the M-functions of Stapp. The covariant formalism 
of the M-functions developed by Stapp [El HE] is based on the 
matrices cr M = (1, (?) and cr^ = (1, — a), where a are the standard 
Pauli matrices. With an each kinematical momentum P M of the 
reaction two by two matrices P = P^o^ and P = P^a^ are as- 
sociated. The products PiPj of these matrices are the elements 
from which the M-functions are built. The matrices Vi and Vi 
associated with the four-velocity of the i-th particle, serve as the 
metric tensors of the particle when performing the contraction 
over this index (traces, successive processes). 

The outline of the paper is as follows. The Sections 1-3 cover 
the three main mechanisms involved in the deuteron break-up: 
Impulse Approximation, NN double-scattering and A-excitation. 
The polarization observables obtained with the covariant deuteron 
density matrix are described in the Section 4. The discussion of 
the calculation results and summary is given in the Section 5. 
The introduction to the covariant formalism of the M-functions 
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is given in the Appendix. 



1 Impulse Approximation 



Pd, &p&n 




Po = ("i,0), r ,a 



P3, T3, 03 (undetected) 

Pi, Ti, Cl 

(spps) 

P2, r 2 ,a 2 (PS') 



Figure 1: IA mechanism. In Saclay experiment pi was measured by the 
magnetic spectrometer SPES-4 in coincidence with the recoil protons (p 2 ) 
detected by Recoil Spectrometer RS. 

The general expression for the S'-matrix elements of the deuteron 
break-up reaction is 

S = i(2tt) /l 5 /l (p 1 + p 2 + P3 ~ Po ~ Pd)M , 

where M is the amplitude of the reaction and we assign indices 
for the each particle of the reaction d(p,2p)n as 

d + p^p + p + n 
(d) + (0) - (1) + (2) + (3) . 
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To avoid in what follows the numerous repetitions of the sum- 
mation sign E, which is inevitable when considering the spin 
and isotopic spin variables, we will use the tensor-like notation 
with the indices characterizing spin and isotopic spin projec- 
tions. The upper and lower indices relate to the final and initial 
channels, respectively. The same index at the up and down po- 
sitions stands for a summation over this index. Free nucleon is 
characterized by the projections of the spin a, of the isospin r 
and by the momentum p. We will use the representation of the 
deuteron spin states by the symmetric spin-tensor S an<7p so that 
the amplitudes of the reaction with participation of the deuteron 
will have the pair of indices a n a p and will be symmetric with re- 
spect to their interchange. 

The amplitude of the deuteron break-up reaction with the all 
indices looks like 

The Pauli principle requires that this amplitude should be an- 
tisymmetric with respect to the interchange of any two final 
nucleons. For example 

m™:37(pi,P2,P3;po,^) = -M;i T ^:2Z 3 (p^PhPs;Po:Pd) ■ 

To fulfill this requirement we proceed as follows. Having written 
the expression antisymmetric with respect to the interchange of 
for example the first and the second nucleons M* 12 ) 3 , we per- 
form then the cyclic permutation of the nucleons. The resulting 
expression M^ 12 ) 3 + M^ 23 ) 1 + M^ 31 ) 2 possess the required antisym- 
metry. 

Let us start with the expression corresponding to the graph 
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pictured in Fig. [H With the all indices it looks like 

%*SS(ft,Mfli,ft) 2m (m * _ m) D T - Ts; ^(p«.K;p<i) . 

(1) 

where ml — (pd — ps) 2 is the squared mass of the virtual nucleon 
and M N n and D are the nucleon-nucleon amplitude and the 
deuteron vertex, respectively. The above expression corresponds 
only to the first term in the development of the spinor particle 
propagator in the Dirac particle-antiparticle formalism 



+ m 1 



m 2 2W 



u ■ U V ■ V 



W = 



W-m W + m. 

Thus we do not consider in what follows the virtual antinucleon 
states, which requires a knowledge of the antinucleon compo- 
nents of the nucleon-nucleon amplitude and of the deuteron ver- 
tex. The deuteron wave function is a product of the vertex and 
the propagator 

2m v (m v — m) p " 

where e is an antisymmetric two by two tensor in isotopic space. 
The fully antisymmetrized 'one nucleon exchange' (ONE) am- 
plitude is equal to 

M ^r 3 :^(pi,P2,^;Po,^) = 
mNNZ%l?Spi^p^Pd - P3y vT3 K;2(pd - P3,P3-,p d ) + 

M^?^?,(P2,P3;po,Pd " Pi)e T " n ^? B (p d " Pl,PlJPd) + (3) 

where we have taken into account the symmetry properties of 
nucleon-nucleon amplitude in the expression (CD). 
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The IA stands for the deuteron break-up amplitude repre- 
sented by only the term in the eq. (jBJ) with the minimal virtual- 
ly of the intermediate nucleon. It corresponds to the minimal 
momentum of nucleon spectator in the deuteron at rest frame. 
Under the kinematical conditions of the Saclay experiment it 
was the undetected neutron (^3). 

In the next subsections we describe the nucleon-nucleon am- 
plitude and the deuteron wave function. 



1.1 NN amplitude 



Pv, C v 



Pi, n, cti 




Po, To, cr 



P2, r 2 , <J 2 



Figure 2: NN amplitude 
The spin and isospin dependent NN amplitude looks like 

pTl„T 2 _ p T 2p Tl Tl T 2 _|_ p T 2p Tl 

2 NL ^ a v \PuP2,P0,Pv ) H 2 NL l ao a v \PhP2iPO,Pv) , 

where e is the unit two by two operator in the isotopic space 
and the Mo and Mi are the isosinglet and isotriplet parts of the 
NN amplitude, respectively. According to the Pauli principle 



8 



the amplitudes Mq and Mi obey the symmetry relations 

M 0aoa'bl^2;^0,^) = M(£j£ (j?2, Ph PO, Pv) , 

MiZl(PhP2;po,Pv) = -Mi^(p 2 ,Pi;po,P,) • 

The cm. canonical amplitudes, i.e. the iS-matrix elements, of 
the NN scattering are expressed in terms of five complex ampli- 
tudes a, 6, c, d, e [19| 

1 A . 

5 = -[(a+6)+(a— 6)n®n+(c+rf)?fi®m+(c— <i)/0/+e(e®n+n®e)] , 

where 

& f k{ t* kf -\- ki _, ki x kf 
m = -=r — , I = -=r — , n = — =r- , 

/cj j are the initial and final cm. momenta of the NN interacting 
system and n = ft ■ a and so on. The shorthand A <g> B stands 
for A^B?. 

<T V (To 

We have to deal with the NN amplitudes in the laboratory 
frame and it is inconvenient to make the transformations from 
numerous individual NN cm. frames to these frame during 
the calculations. It is why the applying of the M-functions of 
Stapp [TZI HE! i s ver y natural in such calculations. The basic 
M-functions analogous to the S'-matrix basis e <8) e, to <S> to, I C*D 
l,n®n,e®n + fi<g)e are built from the products VjYj of the two 
by two (hermitian) matrices V{ = V^a^ and Vj = V- 1 where 

are the four- velocities of the scattered nucleons completed by 
the four-velocity V of the whole system in the arbitrary frame. 
The possible M- function's basis 6j, i = 1, . . . , 6 is given by the 
eqs. (EBHSI1) of the Appendix. Once fixed, the basis allows repre- 
sent the M-function of the NN scattering as a sum 

M = #1 &i + g 2 b 2 + gzh + g^h + g ^ , (4) 
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where gi are the five complex amplitudes. The relations between 
the M-function's amplitudes gi and the canonical amplitudes 
a, 6, c, <i, e are j^Uj 

c + d 

91 = -— ' 

b + a cos <p — ie sin <p 

93 = ^ ' 5,4 

a sin cp + ie cos <£> 



where cj = e l ° , ^ being the cm. scattering angle, and 

(V,V 2 ) + 1 (V,Vo) + l 
(V,V 2 )-1 {V,V Q )-1 ' 

The values of the amplitudes a, b, c, d, e for given laboratory 
energy and cm. angle were calculated by use of the PSA of 
Arndt et al. |21j. The normalization of the M-functions is such 
that the cm. cross-section is equal to 

da 1 Tr(M V V ®V M+ Vi ® V 2 ) 

dtt NN ~~ (8tt) 2 s 4 

The presence of metric matrices Vi and Vj is the peculiarity of 
the Stapp formalism as was mentioned in the Introduction. 

As to the virtual nucleons they are off-shell and we took it 
into account only kinematically. In particular we calculate the 
basis M-functions b{ by use of the eqs. (l28H3TI) using the 'virtual' 
velocity V v = (pd — P3)/m v , if of course ml = (pd — P3) 2 was 
happened to be positive. When it was not the case we assigned 
to the matrix elements of the corresponding nucleon-nucleon 
amplitude zero values. For the scalar amplitudes gi the only 



c — d 



92 = 



2 ' 

b — a cos cp + ie sin ip 



2 

Uq — CO 



UJq — U 



-1 
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way to obtain the off-shell values is the model calculations and 
they are permanently in progress (see for example |22j). Still 
the situation is far from being satisfactory, especially for the 
energies above several hundreds MeV. Thus only the on-shell 
amplitudes obtained in the PSA are of the practical use. The 
choice of the on-shell kinematic corresponding to the off-shell 
one is ambiguous. One way is to take the PSA solution at the 
s = (p + p v ) 2 = (pi + p 2 ) 2 and t = (pi - p ) 2 = (p 2 - p v ) 2 . It 
means to take gi(s,t,ml) = gi(s,t,m 2 ). Another way is to put 
at first the virtual nucleon on the mass-shell p v = (E v ,p v ) — > 
pi = (Jp 2 , + to 2 , p v ), and then get the PSA solution at the 

s * = (Po+Pl) 2 — s an d t* = (p2 — Pi) 2 , to accept that 
gi(s,t,m 2 ) = ^(s*,t*,TO 2 ). We preferred the first way, since it 
ignores the off-shell mass dependence at all in contrast to the 
non-dynamical recipe for such dependence induced by the second 
way. 



1.2 Deuteron wave function 



Pa, c^3 




Pd, CTp 



Pv, <?: 



v 



Figure 3: DNP-vertex 
The DNN-vertex (see Fig. EJ) relates to the wave function by 
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the eq.(J2J). The antisymmetric tensor 




e r v r 3 



determines the isosinglet nature of the deuteron. The classic 
deuteron wave functions are related to the deuteron at rest 
frame. The most suitable for the transition to the M-function 
expression looks like 

$ s ^ 3 (£) = a(e;X 3 + CO - btftfiegef* + efjef*) , 

where subscript s reminds the ^-matrix origin of the classic wave 
function and the unit vector k is directed along the nucleon 
momentum. The matrices e^, i = 1,2,3 coincide with the Pauli 
matrices, but we use the different notation to distinguish them 
from the matrices a 1 , which have another spinor indices: a\ d . 
The scalar functions a and b are connected with the S- and 
.D-wave functions u and w as following 

b b^8 
u = a- - , w = — . 

The corresponding M-function of the deuteron wave function in 
an arbitrary frame is built from the products V v V d and V^V d of 
the two by two matrices, where V£ v 3 are the four- velocities of 
the deuteron, virtual nucleon and on-shell nucleon, respectively. 
It is equal to 

^ (e + V v V d )%(e + V 3 V d )Z + (e + V v V d )Z(e + V 3 V d )% | 
2y/((V v ,V d ) + l)((V 3 ,V d ) + l) 

12 



b (e ~ Wd)g(e ~ VzV d )Z + (e - VW)£(e - 
V((K,^)-1)((^3,^)-1) 



The following normalization equation holds 

tt{^ v*®v v ) = K;?A pdp andn ^lfv,- a ^v v - av(Jv = « 2 +- 2 



where po is the invariant density matrix of the deuteron defined 
in the Section 4 (see eq.(!T§j)). The functions u and w obtained 
with the Paris j23j or Bonn |24| potentials are approximated by 
series of the poles 

u(m% m 2 v ) = ^%Trm d N s \ 9 1 9 - E 



q 2 + a 2 i q 2 + a 2 
w{m% m 2 v ) = V87rm d N D ( * 2 - E 2 1' 2 ) • ( 7 ) 

The values Ns,d are the normalization constants (iVj ~ 0.16 GeV) 
and the dependence on the and m 2 , goes through the 3- 
momentum of the nucleons in the deuteron at rest frame 



q 



m d + m) 2 — m 2 ,) \{m d — m) 2 — m 2 



2m d 



The binding energy e relates to the a as a 2 = em. The sum 
rules Ei = E* d» = 1 should be fulfilled. The dimension of the 
wave functions u and w is GeV^ 1 and they are normalized as 
follows 

J dq(u 2 + w 2 ) = {2ir) 3 2m d . 
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Pd, cr p a n p s ,r s ,a s P3,t 3 ,(T 3 




Po,r ,a P2,r 2 ,cr 1 



Figure 4: Double-scattering mechanism 

2 NN double-scattering 

Let us start with the graph shown in Fig. HI referred to in what 
follows as the DS graph 2(31). Two other graphs with the cycli- 
cally permuted nucleons will be referred to as the DS graphs 
1(23) and 3(12). For the DS graph 2(31) the amplitude with 
the 'spectator' on the mass-shell looks like 

A/r 2(3i) r i T 2 r 3^i^2CT3 . dp s 

M DS r Q ]aoap JPl,P2,P3;PO,P d ) = {27T) 3 2Es 

m™3£(p3,pi;p.,^ 



2mf(rrif — m + ie) 



where p 31 = p 3 + pi = (E 3h p 3 i), E s = ^Jm 2 + pj? and m) = 
(P3i — Ps) 2 - The simplifications are necessary to calculate of 
the integral (JHj) since it requires a knowledge of the off-shell 
nucleon-nucleon amplitudes. The most simple method is to take 
the nucleon-nucleon amplitudes out of the integral sign. Rea- 
sonable to do it at some momentum p s placed on the singular 
surface of the integral corresponding to the mass-shell of the vir- 
tual nucleon /. With the z-axis directed along the momentum 
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P31 the equation of this surface looks like 

r,X n Z _ |P3i| 

(Ps\1 , ( Ps \2 , (Vs 9 \2 



+ (-r + e- 1 —^) 1 = i , 0) 



931 931 931 



W31 



where W31 = v /sii = is the invariant mass of the 31-pair of 
the nucleons and 531 = VS31 ~~ 4m 2 /2 is their cm. momentum. 
At the end of this Section and in the Section 5 presenting the 
calculation results we return to the problem of the choice of this 
Fermi momentum. 

When the nucleon pair at the final state has the invariant 
mass near the nucleon-nucleon threshold it is necessary and 
possible to take into account the off-shell behavior of the corre- 
sponding NN amplitude in the eq.flHD to fulfill the closure sum 
rules in the 1 S% state (see an argumentation in [23 [2SI I2Zj ) . Near 
the threshold it is possible to describe this off-shell behavior by 
the simple form-factor M ^ = M on f(s^i,m 2 ), which is related 
to the 1 Sz wave function of the deuteron (jZj) as follows 

/(*3i, m)) = l- (q 2 + a 2 ) £ -^- 2 , (10) 

where the dependence on the virtual mass and energy goes from 
the cm. momentum 



[(W31 + m) 2 - m 2 f ] [(W31 - m) 2 - m 2 f 

Q = — . 

We used the form-factor (ITOj) up to the energy 200 MeV and 
above this energy we were replacing it by the unity. Note that 
at the low energies the corresponding DS graph is refereed to as 
the final state interaction (FSI) graph. 

Taking the NN amplitudes out of the integral sign, but keep- 
ing inside of it the off-shell form-factor, deuteron wave function 
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and the propagator, we derive the following approximation of 
the DS amplitude from the eq.( 



-^^M^^ s (p 3 ,Pi;p^3i - rf)M^2fe,P3i - P° s ;Po,P d - p° s ) 

F^(s 31: t 2 ) : (11) 

where 

F a v a s( , \ = / d Ps ^ n ( Pd -p s ,p s ] Pd )f(s 31 , (P31-P 
apUn 31, 2j - y 2771/(771/ - m + fe) 

and ^2 = (P2—P0) 2 = (Pd—P3i) 2 - The deuteron wave function $ is 
defined by the eq. (EI) • To preserve the covariance it is necessary 
to approximate the velocities of the virtual nucleon (v ) and the 
nucleon 'spectator' (s) in the latter integral by the corresponding 
values in the NN amplitudes taken out of the integral sign. Then 
the expression for the F becomes 

™ ff , + , F a (e + Wd)%(e + V s V d )Z + (e + Wg;(e + 



2 2 v /((K,v d ) + i)((K,y,) + i) 

F & (e - V v V d )%(e - V s V d )Z + (e - iW)%(e - TO)g 



where 



2V((K,^)-i)((y s ,y d )-i) 



T/ Pd-p° s T/ T/ ^ 



{Pd-P° s ) 2 m m d 

are the four-velocities of the corresponding particles. The com- 
plex functions F a ^ are equal to (see eqs.(E|) 



T? 17 F ^ Ci _ 3Fl 
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where 



dp s u(m 2 dl (p d - p s ) 2 )f(s n , (p 3 i - ps) 2 ) 



{2nf2E s 



2mf(mf — m + ie) 



and F w is determined analogously by the w-wave function. For 
the pole expansion of the functions u and w (EJ) a good ana- 
lytical approximation of these integrals have been obtained by 
J.M.Laget " 



TP ~ 
1 u — 



N s 



-{ 



+4 



quV^2wmd 

P+ P- ( P+ P-\ 

arctan h arctan 2^ Q arc tan h arctan — 

a a i \ a{ on) 

931 , 931 



52 Ci arctan 



In 



2 (a + Oi) 
p 2 _ + a 2 
p\ + a 2 



52 c i c j arctan 
y 2{a i + a j ) 

p*L + a\ 



-{12) 



52 °i In 



P+ + OL\ 



where 



931 = 



[{W 31 + m d ) 2 - t 2 ][{W 31 



m d/ 



2m d 



is the momentum of 31-pair in the deuteron at rest frame and 

S31 + m 2 d - t 2 V 7 ^! 



_ 931 

^ 2 



4m 2 



The replacements in the eq.f ll2|) of the normalization constant 
Ns by the and of the parameters q by the d{ entail the 
expression for the F w . We have tried in the calculations the F w 
thus obtained in addition to the F U: but in the final calculations 
it was still omitted. The reason for this lays in the closure sum 
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rule, which requires in case of the F w ^ to take into account 
the off-shell behavior of the FSI nucleon-nucleon amplitude also 
in the 1 i?3-state, rather than in the ^-state only as we did. 

The real part of the factor F in the eg. (1121) has two terms. 
The first one is due to the off-shell states contribution inside 
the triangle loop, which is proved to be canceled in the eikonal 
regime [T3j. The second term in the real part ( 1121) originates 
from the threshold form factor ( 1101) and we were dropping it for 
the energies of the finally interacting nucleons above 200 MeV. 

Now let us return to the problem of the choice of the Fermi 
momentum p° s in the eq.( llll) . We have considered two candi- 
dates. The first one is advocated as follows. The deuteron wave 
function in the integral (JED drops quickly with the decreasing of 
the virtual mass ml = (pd — p s ) 2 suppressing by this the contri- 
bution of p s corresponding to the small virtual masses m v . It 
is reasonable then to test the p Q s giving the maximum value of 
m v . In the cm. frame of the 31-pair the momenta placed on the 
ellipsoidal surface (JH1) look like (^-,g3in), where n is an unite 
vector. Then 

m l = (Pd ~ P° s ) 2 = mj + m 2 - E d W 13 + 2p d q 31 x , 

where (Ed, p d ) is the deuteron four-momentum in this frame and 
x is the cosine of the angle between n and p^. Therefore the 
maximum m 2 is achieved at x = 1 and the momentum of the 
nucleon s should be directed along the deuteron momentum in 
the 31-pair cm. frame. We will denote the momentum thus 
obtained as p™ ax . The argumentation for the second candidate 
is less evident and consists of the following. Since the states 
with small virtual masses are still contribute to the integral (JED 
we should not use the Fermi momentum corresponding to the 
maximal m v but try a momentum corresponding to some smaller 
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value of the virtual mass m v . We have chosen for this purpose 
the momentum placed at the 'top' of the ellipsoidal surface (El), 
which corresponds to the maximum velocity of the nucleon s 
in the laboratory frame. We will call this momentum as the 
optimal one and will denote it as p° pt . With the z-axis directed 
along the momentum P31 its components are 

P r = (£,o,o,^ +to g), 

whereas the momentum p™ ax has in this frame the components 

Ipl -&L3 
p max = ^ q ^ n ^ q ^ 13_ + 13 n J ? 

Z KK13 

where n is the unite vector directed along the deuteron momen- 
tum in cm. frame of the 31-pair. 

The comparison with the experiment of the calculations per- 
formed using the both Fermi momenta is presented in the Sec- 
tion 5 and definitely testifies in favor of the optimal Fermi mo- 
mentum p° pt . 



3 A -excitation contribution 

The contribution corresponding to the graph shown in Fig. is 
equal to 

M A r ;<Jo(Jp(Jn (Pl,P2,P3;P0,Pd) = 
* 2 2 ' <1=P0-P2i ( 13 ) 
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Ps > T~s i 



p 3 , r 3 , o" 3 (undetected) 

Pi, n , a, (SPES) 
P2, r 2 , a 2 (RS) 



Figure 5: A excitation mechanism 

where the fi is the pion mass. The 7rNN-vertex in the eq. (1131) is 
equal to 

where £ is the isovector index of the pion and 

e +i = (0 o\ 1= /o i\ e „ = Wi o 

U o J ' V j ' 72^0-1 

The spatial part of 7rNN-vertex with the virtual pion looks in 
Stapp formalism like (see the eq.( 12BI) of the Appendix) 

G%(p 2 ,q;Po) = f(q 2 )gMe ~ V 2 V )% , 

where g w ~ 13.6 is the 7rNN coupling constant, and f(q 2 ) is the 
pion form factor, for which we took the monopole representation 

with the cut-off A = 1.0 GeV. 

The spin and isospin dependent 7td — >NN amplitude in the 
eq.f THl looks like 

Mr 3 ;^:bi^3;g,w) = er 3 M^ n (pi,p^q,Pd) , 
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where 

e+1 = ( o o ) ' e - 1 = ( -1 ) • e ° = "72 ( 1 ) ' 

The amplitudes of the physical reactions are equal then to 

M 

We have used for the deriving of the M-function of the ird ^NN 
reaction the ^-matrix elements obtained by J.M.Laget [29j. In 
this approach the one-loop box diagrams with the NA as the in- 
termediate state had been calculated numerically, including the 
7rN intermediate scattering in the S, P and D waves parametrized 
by their phase shifts. However, to avoid the double counting 
with the DS term we have excluded the P\\ -wave (the nucleon 
pole in the 7rN amplitude is a part of the DS term). Finally, 
the p -exchange was considered in the ird ^NN amplitude. The 
details are given in the Appendix of the paper |29j. These cal- 
culations have been performed in the deuteron at rest frame. 
So having at our disposal the ^-matrix elements 

where the subscript (d) means that the momenta are consid- 
ered in the deuteron at rest frame, we have to obtain the M- 
function in an arbitrary frame. The only way in this case is the 
straightforward applying of the common recipe prescribed by 
the eq.(J2BD of the Appendix. At first we derive the M-function 
in the deuteron at rest system: 

M^ n (pi {d) ,P3(dy,q(d),m d ) = v la ^v 3 l 3 S^ n {p 1{d) ,p 3{d )]q {d) ,m d ) . 

(15) 



21 



The boost matrix v% of the nucleon from the rest to the four- 
velocity (Vo, k^/V^ 2 — 1), where k is the unit vector along the mo- 
mentum, is determined by the eq. (1241) . Having obtained the M- 
matrix in the deuteron at rest frame we applied to the eq.tlTHl) the 
boost transformation corresponding to the deuteron moving 
along the z-axis with the four-velocity (Vdo, 0, 0, \jV d \ — 1) so 
that the final M-function looked like (see the eq. (121)1)) 



4 Deuteron density matrix and observables 

Let us remind at first how like looks the density matrix of a 
spinor particle in the Stapp formalism 



p al = l -(V» + S*)of , (16) 

where the a, 6 = 1,2 are the spinor indices, is the four- 
velocity and is the four-vector of polarization of the particle, 
which is orthogonal to V: (V, S) = 0. In the index-less form it 
looks like 

P=\(V» + S»)d,= l -{V+~S). 
The analog of the expression (JTSj) for the vector particle is 

p<J p a p (J n d n _ pfJ,V gOpOp gOnOn (17) 

where 



1 

H 12 
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The S and T are the vector and tensor polarizations with the 
properties 



(s, v) = o , t» v = t v » , VuT^ = o,t; = o. 



In the index-less form eq.(IT7j) could be written as 

P 



(T, 



If the quantization axis is directed along the y-axis then in the 
cm. frame we have 

S =(0,0,^,0), 



/ 



T = 



Pyy 
2^6 



1 





\ 

-2 
I j 



where p y = n + — n_ and p yy = n+ + n_ — 2no- The statistical 
weights n + ,n_,no of the states with the definite projections of 
the spin on the quantization axis are normalized so that n + + 
no + n_ = 1. Then Tr(p 2 ) = n 2 + + Uq + n 2 _. If the particle moves 
along the z-axis with the velocity V = (Vq, 0, 0, V z ), then S does 
not change, but T becomes 



( V 2 

z 



T = 



Pyy 
2^ 





1 
0-2 

\V V Z 



v v z \ 




* J 



V 2 



Pyy 



2V6 



tu,i/ 



We can represent then the p as 



1 



Pd = Po + ^PyPy + ^PyyPyy 



(18) 
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where 

p = -V^V-—a^a^ (19) 
is the unpolarized density matrix and 

p y = -\(°y ®V + V ®a y ) , 

Pyy = ®° V 

are the vector and the tensor polarization matrices, respectively. 
If the (target) proton is unpolarized, the initial density matrix 
of the dp system is equal to 

Pi = ^® Pd 

and the vector and tensor analyzing powers are equal to 
. Tr(M V <8> p y Mt V 1 ® V 2 ® V 3 ) 

^Ay = ■ , (21)) 

. Tr(M V Q p yy Mt Vi <g> V 2 ® V 3 ) 
a A yy = o ' \ 2L ) 

where 

2 

is the unpolarized cross section. The polarization of the fast 
outgoing nucleon in the y direction is determined by the equation 

aP ly = Tr(M p % a y ® V 2 ® V 3 ) , 

where the polarized cross section is equal to 

3 1 

a = Tr(M p t M f V 1 ®V 2 ® V3) = <7 (1 + -PyA, + -^yjAy) . 
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Tr(M V S 


$Py 


Aft 


Vi $ 


$V 2 $ 


^3) 


Tr(M V & 




2 

Mt 






5 V3) 


Tr(M V ( 


3> ^0 


2 

Mt 









The last equation follows from the eqs. (II 802011211) . Introducing 
the polarization Pq of this nucleon for the unpolarized deuteron 
and the depolarization parameters D v as follows 

_ Tr(M V ®po Mt (tV®V 2 ®V 3 ) 

C(K0 = 



2 

Tr(M V ®p y Aft a y ®V 2 ®V 3 ) 

(y^D v = 

we can write polarization P\ y for the vector polarized deuteron 
beam in the following form 

Ply = P ° + } PyDv . (22) 

1 + T^PyAy 



5 Calculation results and discussion 

Let us start with the main parameters of the detecting system 
of the performed experiment. Polarized deuteron beam of the 
accelerator Saturne was incident on a 4 cm thick liquid hydro- 
gen target. Fast protons from the reaction (pi) were selected 
at the scattering angle 6i = 18.0 ± 0.5° by the magnetic spec- 
trometer SPES-4 in the coincidence with the recoil protons (p 2 ) 
detected at the angle 62 = 57.0° (±4.25° in the scattering plane 
and ±8° in the vertical plane) with a mosaic of E and AE scin- 
tillation counters. Two multiwire proportional chambers at 1.5 
and 3 m from the target provided the recoil proton track posi- 
tion. Six central momentum settings for the proton detected in 
SPES-4 O10 = 1.6, 1.7, 1.8, 1.9, 2.0, 2.05 GeV/c) were studied. 
The polarimeter POMME located behind the final focal plane of 



25 



-2 
-2.5 



i_ (ONE) 




Tkl 7. tdJ 


Tk31,td2 
~ ... I ... I 


Tk23, td' 
. . . 1 . . . 1 . . . 


1 . . . 1 . . . 


0.2 0.4 


0.6 0.8 


1 1 .2 




TDI VS TKIJ 


Tk,j. GeV 


f (DS) 

^ Tk31,_t02 

" , , , 1 , , , 1 




V 


Tk23, -t01 

. . . I . . . I . . . 


Tk12, -103 



0.2 0.4 0.6 0.8 1.2 

Tk,:. GeV 

TOI VS TKL 



Figure 6: The (Tij,t di ) and (Tij,t 0i ) distributions for p 10 = 1.6 GeV/c and 
HES. 

SPES-4 was used to measure the polarization of the fast protons. 
At the chosen Oi and for the momentum pi of the fast proton, 
two kinematic solutions are possible for the recoil proton energy 
T2, which according to our definition are the high energy (HES) 
and the low energy (LES) solutions. The neutron spectator mo- 
menta in the deuteron rest frame, denoted as q, ranged from 30 
to 440 MeV/c. 

Thus the calculation procedure should include scanning over 
the acceptance of the detecting system. By scanning over the 
allowed phase space of the detectors, we intended to obtain a re- 
alistic calculation result, which can be directly compared to the 
data. This scanning procedure inevitably invokes the necessity 
for a Monte-Carlo type event generation, because the number 
of phase unit volumes amounts to a huge number when the unit 
volume is defined by the resolution of the measurement. 
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We would like to show now, taking the p\o = 1.6 GeV/c 
and HES as an example, some kinematic distributions charac- 
terizing different mechanisms of the deuteron break-up in this 
experiment. In Fig. El the two-dimensional distributions of the 
kinetic energies = (sij — Am 2 ) /2m of the i,j pairs of the 
final nucleons with two momenta transfer, = (p^ — Pi) 2 and 
toi = (Po~Pi) 2 , are shown. Note that tdi are equal to the squared 
masses of virtual nucleon in the nucleon pole graphs. It is seen 
that of three final nucleons considered as the spectator in these 
pole graphs, the neutron (^3) provides the maximal masses of 
the exchanged nucleon. The pole graph with P2 as the spec- 
tator gives the negative squared masses and we have excluded 
this graph from the calculation. The distributions of (T^-, toi) 
characterize the graphs of the DS mechanism. T31 ranges from 
30 to 200 Mev and therefore the DS graph 2(31) is the typical 
FSI graph. Two other pairs of the final nucleons and especially 
two protons pi and P2 have much higher relative energies. How- 
ever applying of the Glauber approach when calculating the DS 
graphs with these pairs being rescattered (121 Ej is not justi- 
fied because the momenta transfer £93 and toi for these graphs 
are very high. Thus we have applied for calculation of the DS 
graphs the expression ( 1111 ) and two others obtained from it by 
the cyclic permutation of the final nucleons. 

We would like now to return to the problem mentioned in 
the Section 2. It is the choice of the Fermi momentum p® in 
the eq.flTTj). The dotted line in Fig. presents the polarization 
observables calculated with the p™ ax in the DS graph 1(23) and 
with the p° pt in the DS graph 3(12) (the sensitivity of the con- 
tribution of the FSI to the choice of Fermi momentum is weak 
and we fixed it to be p™ ax ) . The dashed line corresponds to the 
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Figure 7: Polarization observables for the different choices of the Fermi mo- 
mentum. Dotted line presents the calculations with p™ ax in the DS graph 
1 (23) and p° pt in the DS graph 3(12). Dashed line corresponds to the choice 
of p° pt for DS graph 1(23) and p™ ax for the DS graph 3(12). Solid line is 
obtained with the p° pt in the both these graphs. The dashed-dotted line 
presents the only FSI prediction (DS graph 2(31)). 
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choice of the p° pt for the DS graph 1(23) and of the p™ ax for the 
DS graph 3(12). The solid line is obtained with the p° pt in the 
both these graphs. The improvement in the description of the 
data in the latter case is evident. It is important that neglecting 
of these graphs results in the significant worsening of the data 
description. The dashed-dotted line presents the only FSI pre- 
diction and we see that it does not provide the suppression of 
the tensor analyzing power at the high q observed in the experi- 
ment. So we can summarize that besides the FSI graph the both 
1(23) and 3(12) DS graphs are required for data description and 
they should be calculated with p opt as the Fermi momentum. It 
is the very choice of the DS model, which we have applied in 
further calculations. 

The measured tensor analyzing powers A yy and the calcu- 
lation results with Bonn deuteron wave function are shown in 
Figs. El and El The main feature of the experimental data is the 
strong deviation from the IA (dashed line) for q > 0.2 GeV/c. 
The full theory including the A-excitation graphs (full line) gives 
the satisfactory explanation of this deviation, the DS graphs 
(dashed-dotted line) playing the decisive role. One could expect 
the more noticeable A-effects taking in mind that the laboratory 
kinetic energies of the 23-pair of the final nucleons range from 
0.7 to 0.9 GeV (see Fig. ED, which is exactly the region of the 
A dominance in the ird -^NN amplitude. Yet the momentum 
transfer £qi from the target proton to the fast proton detected by 
the SPES-4 is very high and the pion form factor ( TH1) reduces 
the common contribution of the graph shown in Fig. [5] ( with, 
of course, correspondingly interchanged final nucleons). 

The vector analyzing power A y results, shown in Figs ITU and 
ITT1 exhibit the similar tendency: significant correction to the IA 



29 




Figure 8: The tensor analyzing power A yy for the HES. The experimental 
points are presented for the different values of the central momentum detected 
in the magnetic spectrometer and as a function of the outgoing neutron 
momentum in the deuteron at rest frame. The dashed-doted line is the IA , 
the dashed line has in addition the DS contribution, and the continuous line 
is the full calculation including in addition the virtual A. 
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Figure 9: The tensor analyzing power A yy for the LES. Same notations as in 
Fig. El 
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Figure 10: The vector analyzing power A y for the HES. Same notations as 
in Fig. E 
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Figure 11: The vector analyzing power A y for the LES. Same notations as 
in Fig. E 
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above 0.25 GeV/c by the full calculations. The surprisingly good 
coincidence, taking in mind the very small statistical errors of 
the measured values, of full calculations with the experimental 
points is achieved at the all SPES-4 settings for the both high 
and low energy solutions. 

For the polarization of the fast proton the data are averaged 
over the SPES-4 settings to obtain more statistics. However the 
calculations were performed for each setting. The experimen- 
tal data and the calculations are presented in the terms of the 
proton polarization Pq for the unpolarized deuteron and the de- 
polarization parameters D v (see the eq. (l22l) ) and are shown in 
Figs fHflTHl and respectively. For HES a good description is 
obtained for the P$ with the full diagram calculations, whereas 
it is not the case for LES. For the D v , the IA is even closer 
to the data than the full calculations. However, considering the 
large error bars for the D V: one can say that there is not decisive 
discrimination between the both calculations. 

The above results were obtained with the Bonn deuteron 
wave function |24j. Returning to the initial idea of this ex- 
periment to discriminate between the different deuteron wave 
functions we have tried also widely used Paris deuteron wave 
function [23j . In Fig. [TBI the results of the calculations for the 
both wave functions at 1.8 and 2.0 GeV/c settings of the SPES-4 
are shown. It is seen that the dependence on the deuteron wave 
function of the full calculations is weak. 

To conclude we would like to remind that the main task of this 
experiment was the investigation of the deuteron structure in 
the high internal momentum region. The polarization data were 
thought to bring the information on the S- and .D-component of 
the deuteron wave function, which would be possible only if the 
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Figure 12: The forward proton polarization Pq for the HES. The notation of 
the curves are the same as in Fig. |H1 The calculations are done consistently 
as for A y and A yy for each setting of the spectrometer while the data are 
summed. 
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ure 13: The forward proton polarization P f° r the LES. The notation of 
curves are the same as in Fig. |HJ 



36 





0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

q (GeV/c) q (GeV/c) 



0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

q (GeV/c) q (GeV/c) 



G 1 




2 GeV/c 


° 1 




2.05 GeV/c 


0.5 




t * 

;:>!. 


0.5 




( ' 















0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

q (GeV/c) q (GeV/c) 



0.1 0.2 0.3 0.4 0.5 

q (GeV/c) 



0.1 0.2 0.3 0.4 0.5 

q (GeV/c) 




0.1 0.2 0.3 0.4 0.5 

q (GeV/c) 



0.1 0.2 0.3 0.4 0.5 

q (GeV/c) 



Figure 14: The depolarization parameter D v for the HES (right picture) and 
LES (left picture) . Same notations as in Fig. [S| 
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Figure 15: Comparison of the Bonn(full line) and Paris(dashed line) deuteron 
wave functions predictions 
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IA would dominate in the reaction mechanism. The underlying 
idea of this experiment as well as the analogous ones is that the 
deviations from the IA, corrected by other presumingly small 
contributing conventional mechanisms, could be interpreted as 
the revealing of the quark degrees of freedom. In fact, the po- 
larization parameters, measured at first time in this exclusive 
experiment with a large accuracy and in a large amount, deviate 
strongly from the IA at the high internal momenta. However, 
the estimations of other traditional reaction mechanisms, the 
nucleon-nucleon double-scattering being the most significant, 
were tried and have shown a qualitative agreement to the re- 
ferred experimental data. This agreement between the data and 
the theory is thought to imply that the possible new degrees of 
freedom in the deuteron structure contribute much smaller than 
the hadronic degrees. Still definite perfection of the model is 
necessary for the adequate description of the polarized deuteron 
break-up data. Investigation of the short range deuteron wave 
function, if it is at all possible in the deuteron break-up ex- 
periments with the hadron probes, will be, by our opinion, the 
further step only after all discrepancies in theoretical description 
of polarization observables are overcome. 

Acknowledgements 

The author thanks S.Belostotski, A.Boudard, J.-M.Laget and 
V.Nikulin for fruitful discussions and cooperation during per- 
forming this work. He is also very grateful to A.Boudard for 
the warm hospitality at Service de Physique Nucleaire, Centre 
d'Etudes de Saclay, where the main part of this work was done. 



39 



Appendix. M-functions 



Let us consider the process with the rif and rii spinor particles 
in the final and initial channels, respectively. We will denote the 
corresponding amplitude with the indices as 



M^ -' T "% ...... 



Let Vi be the boosts of i-th spinor particle moving with the 
velocity V{. Then the connection between the ^-matrix and M- 
function looks like 

S = <g> . . . <g> v~j M Vi . <g> . . . <g> v ni . (23) 



The boost matrix v% from the rest to the 4- velocity (Vb, Vq — 1), 
where k is the unit vector along the momentum, is equal to 

vi =( Vo + V3 Vl - W2 ) , (24) 

\Vi+ IV 2 Vq-V-s J 

where 

^ = -^(v / H+T,k v / TW) • 

With respect to Lorentz transformations, represented by the 
unimodular matrices z 



~a _ 

Z b — 



zo + iz 5 izi + z 2 

iZi - Z 2 Zq- iz 3 



where 20,1,2,3 are the complex numbers such that ZQ+zf+z^+z^ = 
1, the M-functions behave as the operators in the direct product 
of spinor spaces 



z® ...®zMz 1 <g> ...(g) z 1 . (25) 
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The simplicity of the transformation ( 1251 ) is the main advan- 
tage of the M-functions. It allows to build the M-functions from 
products VjVj of two by two matrices V{ — Vfd^ and Vj = Vf<j^ 
where V- 1 are the four- velocities of the scattered spinor particles 
completed by the four- velocity V of the whole system in the arbi- 
trary frame. The matrices = (1, (?) and = (1, — a), where 
a are the standard Pauli matrices, have the different position of 
spinor indices specified in the following way 



° ^ °ab J ° ^ ° 



The simplest M-functions are easy to derive using the Dirac 
bispinors in the Weyl representation in the M-function form 



75 = 



u% = \M yp\ , u b a = V^(e b a v b ~ a ) . 

Here a is the bispinor index, Vab = V^nb j = V fl a^ ba and 
is the four- velocity of the particle. Let us consider for example 
the 7rNN vertex. In the Weyl representation the 75- matrix looks 
like 

' e 
-e 

and the amplitude of the transition between real nucleons with 
the production of the pion is equal to 

M?(p/, q; Pi ) = <?X75?*A = 9M4 ~ VfV m ) , 
or in the index-less form 

M b/, 9; Pi) = 9wUfj5Ui = gnm(e - V/Vi) . (26) 

The matrices Vi and Vi, where V/ 1 is the four- velocity of the i- 
th particle, serve as the metric tensors for the spin index of this 
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particle in Stapp formalism when performing the contraction 
over this index (traces, successive processes). So, in contrast 
to the simple expression of the cross section via the ^-matrix 
elements 

a = —Tr(SS^) , 
the same via the M-functions looks like 

o = ^m^.zZ'v^ . . . v:;^M7/'t f vi f - . --Vn 

or in the index-less form 

a = ^-Tr(M V u ® . . . ® V ni M%. ® . . . ® V nf ) . (27) 

At last let us give one of the basis of the M-functions for 
the NN scattering. It is convenient to use such functions bi in 
the development (Jlj), that the cross section is equal to £i \gi\ 2 - 
From the eq.(l2~7j) follows that for this purpose one should build 
the basis orthonormalized with respect to the scalar product 

AyM SpjMiVo ® VM)Vi ® V 2 ) 
(Mi, Mj) = — J - . 

The eqs.(J2EHSH) present such a basis 



= * = 1 -- -6, (28) 



. JMi,| 
where 

Mi = (e - ViVj ® (e - V 2 V ) . 

M 2 = (ViV - VK) ® (W - VVo) • 

M 3 = (e + V X V V ) ® (e + V^o) , 

M 4 = (a lu (e + ViK) (^i^ + VK)) ® (a 20 (e + ^Vb) - (V 2 V + VV ) 

M 5 = (e + ® (a 20 (e + V 2 y ) - (V 2 V + VVo)) , 

M 6 = (a lv (e + V^K) - (ViV + V%)) ® (e + V 2 V ) , 
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with 



(Vt + Vj) 2 . (30) 



and 



Milr = ±[l-(V h V v )}[l-(V 2 ,V )} 
M 2 \\ 2 = 4 [1 + (Vi, V v ) - 2(V h V)(V V , V)} [1 + (V 2 , Vb) - 2(V 2 , V)(V , V)] 
M 3 || 2 = 4[l + (Vi,K)][l + (V 2 ,V )] 

iVi 4l! - p^p 

M 5ll 2 = {[1 - (V, V 2 f] [1 - (V, V ) 2 ] - {(V 2 , V ) - (V, V 2 )(V, V )} 2 } 

M 6ll 2 = ^{[1 - (y,K) 2 ] [1 - (^K) 2 ] - - (^^)(KK)f} • 

The Vi in the above equations are the four-velocities of the i-th 
particle, V being the four- velocity of the whole NN cm. system. 
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